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There  exist  examples  of  systems  in  which  small,  simple 
metallic  cations  are  associated  with  ring-like  molecules 
(c.g.,  aromatic  radical  anions),  with  aggregates  of  molecules 
which  encage  the  ion  (e.g.,  solvation  and  coordination),  or 
special  molecules  which  trap  ions  within  the  molecular  structure 
(crypts).  For  each  of  these  examples,  the  harmonic  oscillations 
of  the  ion  in  the  presence  of  the  molecular  structure  can  be 


JO  con  t  i  lined 


observed.  Moreover,  in  many  instances,  the  stability  of  the 
ionic  motions  in  the  presence  of  the  molecular  structure  is 
important  in  considering  some  form  of  transport  process  which 
is  associated  with  the  ion.  In  this  paper  we  present  an  analysis 
of  the  mechanics  of  the  motion  of  an  ion  under  the  influence  of 
various  molecular  structures.  for  our  analyses,  we  use  a  form 
of  symmetry-adapted  Taylor  series  which  we  recently  developed. 

We  determine  the  positions  of  equilibrium  for  the  ion  in  the 
presence  of  a  continuous  ring,  a  polygon  and  a  polyhedron  of 
bourses  for  the  pair-wise  forces.  We  also  determine  the  states 
of  stability  for  the  positions  of  equilibrium  through  an 
examination  of  the  second  order  terms  in  the  Taylor  series. 

We  show  that  the  general  form  for  the  Taylor  series  for  a 
continuous  ring  is  of  the  same  form  as  that  for  a  discrete  \ 
polygon  .  We  also  show  that  small  numbers  of  terms  in  the 
Taylor  series  can  duplicate  the  exact  potential  for  an  atom 
within  a  tetrahedron  or  octahedron  of  sources  to  a  high  degree 
o f  accuracy , 
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Abs tract 

There  exist  examples  of  systems  in  which  small,  simple  metallic 
cations  are  associated  with  ring-like  molecules  (c.g.,  aromatic 
radical  anions) ,  with  aggregates  of  molecules  which  encage  the  ion 
(e.g.,  solvation  and  coordination),  or  special  molecules  which  trap 
ions  within  the  molecular  structure  (crypts).  For  each  of  these 
examples,  the  harmonic  oscillations  of  the  ion  in  the  presence  of  the 
molecular  structure  can  be  observed.  Moreover,  in  many  instances, 
the  stability  of  the  ionic  motions  in  the  presence  of  the  molecular 
structure  is  important  in  considering  some  form  of  transport  process 
which  is  associated  with  the  ion.  In  this  paper  we  present  an  analysis 
of  the  mechanics  of  the  motion  of  an  ion  under  the  influence  of 
various  molecular  structures.  For  our  analyses,  we  use  a  form  of 
symmetry-adapted  Taylor  series  which  we  recently  developed.  We 
determine  the  positions  of  equilibrium  for  the  ion  in  the  presence 
of  a  continuous  ring,  a  polygon  and  a  polyhedron  of  sources  for  the 
pair-wise  forces.  We  also  determine  the  states  of  stability  for  the 
positions  of  equilibrium  through  an  examination  of  the  second  order 
terms  in  the  Taylor  series.  We  show  that  the  general  form  for  the 
Taylor  series  for  a  continuous  ring  is  of  the  same  form  as  that  for 
a  discrete  polygon.  We  also  show  that  small  numbers  of  terms  in  the 
Taylor  scries  can  duplicate  the  exact  potential  for  an  atom  within 
a  tetrahedron  or  octahedron  of  sources  to  a  high  degree  of  accuracy. 


1.  Introduction 


There  are  many  examples  of  systems  which  involve  the  associa¬ 
tion  of  an  ion  or  atom  with  a  ring-like  molecule,  with  a  cage-like 
molecule,  or  with  an  aggregate  of  molecules.  Cation/radical  anion 
pairs,  in  which  the  radical  anion  is  an  aromatic  ring,  constitute 
a  class  of  examples  which  has  been  well  studied.*  In  the  last 
decade  or  so,  considerable  interest  has  grown  around  certain  molecules 

referred  to  generically  as  crypts,  which  are  capable  of  binding 
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small  ions  within  a  fairly  rigid  molecular  cage.  Finally,  the 
well  known  states  of  solvation  or  coordination  of  ions  in  solution 
also  serve  as  examples.^  For  all  of  these  systems  of  aggregation, 
it  is  of  considerable  importance  to  be  able-  to  specify  the  states 
of  vibration  and  the  conditions  for  mechanical  stability  for  an  ion 
or  atom  with  respect  to  the  larger  ring,  crypt,  or  cage  of  solvent 
or  ligand.  A  number  of  transport  and  kinetic  processes  depend  upon 
these  vibrational  states  in  order  to  account  for  part  of  their 
act ivation. 

In  this  paper,  therefore,  we  examine  the  vibrations  and 

mechanical  stabilities  of  an  ion  or  atom  which  is  associated  with  a 

ring,  crypt,  or  cage  of  solvent.  In  order  to  carry  out  this 

examination  of  the  vibrations  and  mechanical  stabilities,  we  use 

a  form  of  symmetry- adapted  Taylor  series  which  we  recently 
4-6 

developed.  An  important  facet  of  the  analysis  to  follow  builds 

upon  the  fact  that  polygonal  rings  and  polyhedral  cages  are  closely 
related  from  a  mechanical  point  of  view. 
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In  the  next  section  we  summarize  the  formulae  which  are  needed 
to  carry  out  the  subsequent  analyses.  Following  that  summary,  we 
present  the  results  of  our  examinations  of  the  ring  and  cage  systems. 
We  examine  the  use  of  a  variety  of  potential  energy  functions  from 
the  electrostatic  Coulombic  potential  to  the  Buckingham  exp-6 
and  related  potentials. 

We  have  already  shown  that  the  form  of  analysis  summarized 
in  the  next  section  is  very  useful  for  the  elaboration  of  the  nature 
of  forces  which  operate  in  solvated  ionic  systems. ^  These 
methods,  we  believe,  also  can  be  useful  for  many  other  applications 
as  well. 

2.  The  Symmetry-adapted  Taylor  Series 

Flsewhere  we  have  shown ^  that  for  a  general  class  of  functions 
which  consist  of  separable  and  purely  angular  and  radial  parts  of  the 
form 


(HR)  =  YXv(R)F(R),  (2.1) 

in  which  Y  x  (  R )  is  the  spherical  harmonic  function,  the  Taylor 
series  is 


C(R*r)  -  (4.)3/2  l  £  l  (•i)L*"AntYLM(R)Ytn,(i) 

~  n=0  n>  L,M,£,m  n  L 


f 


2  F,+ 1 


p2£+l )  (2X  +  1)  J 


1/2 


(L£00 | AO) (L£Mm| Xp) Inl (R) . 


(2.2) 
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The  vector  r  is  a  displacement  from  an  end-point  of  R.  In  eq  (2.2) 
the  quantity  A  ^  is  given  by^ 


A  „  =  0  for  l  >  n  and  n  -  £  odd 
n  £ 


(2£+l)n! (n-£+l) ! ! 

In'-  H+DTTn'+r+TJ  n 


for  £  <  n  and  £ 


n  even 


(2.3) 


1 2 

The  quantity  (L£Mm|Au)  is  a  Clebsch-Gordan  coefficient  and  I  j  (R) 
is  defined  by 


r°° 


TnL^  = 


(277) 


dk  kn+2  f (k) jt (kR) 


(2.4) 
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in  which  j^(kR)  is  the  spherical  Bessel  function  of  the  first  kind. 
The  function  f(k)  defines  the  radial  Fourier  transform  of 
F(R) :5 


f  (k) 


4  77  i  ^ 


°n 

dR  R2F(R)jx 

0 


(kR). 


The  Taylor  series  for  a  scalar  function. 


(2.5) 


L 


G0(R)  -  Yoo(R)  [/^F(R)1 
=  F(R), 


(2.6) 


assumes  a  particularly  simple  form: 

(io(r  +  R)  =  /47  l  (rn/n ! )  [(-i)l+nA  P  (R-?)I  (R)  (2.7) 

n=0  £  n  nii 

i  2  *  * 

in  which  P^(x)  is  the  Legendre  polynomial  of  order  £  and  R«r  is 


J 


the  scalar  product  of  the  unit  vectors  R  and  r. 

We  have  also  shown^  that  the  Taylor  series  (2.2)  for  a 
scalar  function  can  he  expressed  as 


F(R+r) 


<.x»  /\  /\  2. 

I  (rn/n!)£.\  P  (R-r)  l 
n=0  «  nv  q=0 


(~l)q(g>q')  ! 

( S  - <1 1  !  ( Iqi  !  ! 


n-q 

K 


+ 


d' 

dR 


(d/dR) n  _q " 1 H (R)  ( 


This  simple  form  retains  the  angular  dependencies  of  eq  (2.2) 
and  displays  differentiations  with  respect  only  to  the  radial 
quant i t i es . 

We  now  list  several  explicit  expansions  of  potential  energy 

functions  which  we  will  use  in  the  next  section. 

The  expansion  for  the  Coulomb  potential  is  identical  to  the 

1  3 

t  ami  1 Lar  Laplace  expansion 


TFrT 


v  (-r) 


i  R 


f+T 


VR 


■r) 


For  the  Morse  potential, 


<t>M(R)  =  Dexp  [a  (Ro -R)  Hexp  [a  (R0 -R)  ]  -  2}  ( 

where  D  is  the  dissociation  energy  and  a  is  a  parameter,  we  have 
f  ound^ 

*M(R*r)  =  n  l  (-ur)n  X-  An)lP)l(R.O{2ne2oRo[2aRkp_1(2aR) 

n  ,  Z 

-  fn-t)kp(2aR)]-  2cu  R°  [a  Rk  f  a  R)  -  (n-ll)k  («R)]  } 


( 
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in  which  the  functions  k^(x)  are  the  modified  spherical  Bessel 
functions  of  the  third  kind:^ 


kn(x)  =  (-x)n(x'1d/dx)n 


(2.12) 


For  an  arbitrary  inverse  power  of  R,  R 


-q 

9 


we  have  also  found 


5 


R+r 


i  l 

(q- 2) ! Rq  n,£,s 


WR*r) 


U+s) ! (q+n-s-2) ! f 

n!Tr-StTT2T5'r  c’r/RJ 


(2.13) 


14 

This  expression  is  equivalent  to  one  derived  by  Briels  who  used 
functional  expansions. 

The  Buckingham  exp-6  potential 


4>g  (R)  =  A  e'BR  -  C/R6 


(2.14) 


is  simply  a  combination  of  an  expansion  for  the  exponential  and  eq 
(2.13).  The  result  is 


<J>B(R+r) 


=  l  C 

n ,  £ 


-vn  1 

r)  rTT 


An*VR‘r) 


ABn [BRk 


t_1(BR)-(n-£)kjl(BR) 


_L  1  V  Cjjs)  ?  (n+4-s) 

24  Rn+6  («,-s)  !  (2s)  !  !  j 


(2.15) 


It  is  a  straightforward  matter  to  add  the  Coulomb  potential 
to  any  of  the  various  forms  in  order  to  generate  potential  functions 
for  ions  which  interact  with  a  variety  of  molecular  aggregates. 

Crystal  field  calculations  with  the  use  of  the  Coulomb  potential 
eq  (2.9)  of  course  are  well  knownd^  An  essential  feature  of  the 


cyrstal  field  calculations  is  the  exploitation  of  cubic  point 
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group  symmetries  to  sort  out  various  terms  in  the  expansions.  The 
symmot ry - adapted  Taylor  scries  we  have  developed  works  similarly. 

The  examples  we  investigate  in  the  next  section  illustrate  this. 

5 .  App 1  Teat  ions  to  Rings  and  Cages 

i 

We  begin  this  section  with  an  investigation  of  the  system  which 
consists  of  a  continuous  ring  of  matter.  We  develop  a  Taylor  series 
for  the  potential  energy  function  for  such  a  system.  Next,  we 
examine  a  polygon  of  sources  in  a  plane  (an  example  of  which  would  be 
the  benzene  ring).  We  show  that  the  results  obtained  for  the  poly¬ 
gon  bear  a  strikingly  close  similarity  to  those  obtained  for  the  ring, 
finally,  we  examine  systems  of  sources  which  are  distributed  at  the 
vertices  of  a  tetrahedron  and  an  octahedron.  We  develop  these  analyses 
respectively  from  cons iderat ions  of  distributions  of  sources  in 
triangles  and  squares.  Additional  sources  are  placed  at  appropriate 
positions  on  the  axes  which  pass  through  the  centres  of  these 
figures.  The  conditions  for  mechanical  stability  within  these 
systems  still  hear  some  similarity  to  the  conditions  for  stability 
which  we  derive  for  the  ring  and  polygon. 

a_.  Ri n g_s  of  so urce  density 

The  classical  problem  of  the  determination  of  the  form  of 
the  potential  due  to  a  ring  of  charge  density  is  well  known. ^ 

Similar  treatments  for  other  possible  source  matter  densities  seem 

not  to  have  been  given.  We  have  considered  this  problem 

for  the  Yukawa  (or  Debye -IKickel )  potential**  and  we  examine  it 
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further  in  this  section. 

We  consider  the  general  problem  of  the  potential  due  to  an 
incremental  source  dQ  distributed  continuously  as  a  ring  of  matter 
(neutral  or  charged) : 

dQ  =  dV  (3-1} 

The  magnitude  of  R  is  a  constant.  Thus, 

dQ  VR*r)  =  {  d4»R  l  Pjf(cos0R)P;(Cos0  )exp[im(d>R 

J  J  o  m 

*  <?>r)] 


Pi (cos  0R) I  1 1 (R)  =  0.  (3.4) 

From  this  equation  we  identify  two  conditions  for  equilibrium: 

(1) 

Pi(coseR)  =0,  0R  =  tt/2  (3.5) 


The  position  of  the  point  of  equilibrium  is  at  the  origin  of  r,  viz., 
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the  centre  of  the  ring. 

The  second  condition  for  equilibrium  is 

(2) 

In(R)  =  0.  (3.6) 


This  condition  may  in  fact  be  satisfied  by  more  than  one  value  of 
R.  Usually,  as  we  shall  see,  only  one  value  of  R  is  significant; 
the  second  value  may  correspond,  for  example,  to  the  (not  useful) 
case  of  infinite  separation. 

It  is  worth  noting  at  this  point  that  from  eq  (2.8) 


In(R) 


dF 

77 


0 


(3.^) 


as  one  would  expect. 

Stability  is  determined,  for  an  arbitrary  form  of  displace¬ 
ment,  by  evaluating  the  second  order  coefficient  in  the  Taylor 
series.  The  second  order  term  from  (3.3)  is  written  as 


t2 (R+r) 


(R)  ♦ 


(COS0R) P2  (cos0r)  I  2  2  (R) 


(3.8) 


For  the  first  point  of  equilibrium,  (1)  above,  0„  =  it/2,  P2(cos0r  )  = 

e  e 

-1/2,  and  R  =a,  the  radius  of  the  ring,  and 


t 


r“ 

7 


•  2  0 


+  P 


2 (co  s  0  r ) r22 


(Re) 


(3.9) 


For  an  axial  displacement  at  the  centre  of  the  ring,  P2(cos0r) 


=  +1.  Thus, 
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(axial ) 
Z2 


/Jl\ 

' y 


:It(Re) 


*  I 


2  2 


(Re) 


(3.10) 


The  sign  of  I2o  +  1 2  2  is  all-important.  If  the  sign  is  negative, 
the  motion  is  stable. 

For  a  transverse  displacement  about  a  position  of  equilibrium 
at  the  centre  of  the  ring,  P2 (cos6r)  =  -  1/2,  and 


(trans) 

z2 


/Tn 


r 

I 


l20(Re) 


j!  2  2 (Re) 


(3.11) 


The  second  position  of  equilibrium  lies  on  the  axis  which 
passes  through  the  centre  of  the  ring.  At  this  second  point,  R=Re 
where  Rg  is  the  solution  to  1 1 1 (R  )  =  0.  We  find, 

cose  =  ^  [R*-a2]1/2  (3.12) 

e  e 

with  the  Legendre  polynomial  evaluated  as 


P2(cos0£  )  =  1  -  3az/2R2 


Thus , 


(3.13) 


1 2  ^  =  /4tt  Q 

For  an  axial 

. (2 , axial) 

Z2 


j[-  jl 2  0 (Rg)  +  f(l-3a2/2R2)P2(coser)I22(Re) 

displacement,  P2(cos0r)  =  +1,  and 

/4¥  Q  j[-  ^2o(Re)  +  (2/3-a2/R2)I22(Re)  . 


(3.14) 


(3.15) 
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lor  a  transverse  motion,  I’pfcosO  )  =  -  1/2,  and 

t (2 , trans)  =  ^  Q  r 2 f  ll2o(Iy  -  i(l-3a2/2R2)I22(Re)j .  (3.16) 

It  is  difficult  to  make  further  progress  at  this  point  without 
considering  specific  functional  forms.  From  the  expansions  of  the 
various  functions  given  in  the  last  section,  we  are  able  to  consider 
several  cases. 

We  examine  now  the  ring  of  source  for  the  Morse  potential. 

The  equilibrium  condition  I i i (R  )  -  0  specifically  yields 

-  2ctDexp(a(Ro-Re)]  exp  [a  (R0  -  Rg)  ]  - 1 J  =  0.  (3.17) 

We  see  that  Re  =  Ro  or  Rg  =  00  satisfy  this  equation.  The  latter 
value  is  not  useful. 

The  quantity  Ro  in  the  Morse  potential  is  a  separation  which 
is  characteristic  of  the  Morse  interaction  which  operates  between 
two  individual,  spherical  masses.  It  is  definitely  the  case  that 
when  the  mass  is  distributed  on  a  ring,  positions  and  conditions  of 
equilibria  and  stability  need  not  conform  to  the  conditions  which 
apply  to  the  simpler  spherical,  diatomic  case.  In  this  respect, 
with  the  use  of  the  Morse  (and  other)  potentials  in  extended  distribu¬ 
tions,  we  see  arising  conditions  of  interpretation  which  are  necessary 
for  composite  potential  functions  such  as  the  Stockmayer  potential.1^1 
The  Stockmayer  potential  is  the  sum  of  a  Lennard- Jones  component  and 
a  dipolar  interaction.  The  Lennard- Jones  component,  operating  alone, 
has  a  meaningful  minimum  value  which  is  associated  with 
a  point  of  equilibrium  for  a  diatomic  system.  When 
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comb  ined  with  dipolar  contributions,  as  is  the  case  with  the 

Stockmayer  potential,  the  minimum  of  the  composite  function  need 

not  fall  at  the  same  place  as  the  minimum  for  the  Lennard -Jones 

1 7 

component  alone.  Yet,  typically  an  interpretation  is  given  which 
says  simply  that  the  Lennard- Jones  type  of  minimum  in  a  Stockmayer 
potential  gives  the  location  of  the  minimum  if  the  dipolar 
component  were  to  vanish.  We  see  similar  interpretations  applying 
to  the  use  of  diatomic  potentials  for  extended  distributions  of 
mass  and  charge. 

The  second  order  term  for  the  expansion  of  the  Morse  potential 
is  simply 


r2  2 


PM2 


—  -jci2D  exp  [a(Ro  -R)  ] 


2  exp [a (Ro -R) ] (1 -1/aR) 


-  1  +  2/aR  +  P2 (cos9r) P2 (cos 0R) 


4  exp[a(R0-R) ] (l+l/2aR) 


2 (1+1/aR) 


(3.18) 


The  individual  contributions  I2o  and  I22  are  displayed  separately 
as 


T (Morse) 
l20 


— —  a2D  exp[a(Ro-R)] 

/TW 


2 


exp [a (R0 -R) ] (1-1/aR) 


and 


1  +  2/aR 


(3.19) 
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.  (Morse ) 
1  >  1 


a2n  exp  [  a  (.Ro  -  R)  J 
/4u 


2  exp|n(R0-K) | (l+l/2uR) 


-  1 


1/«R 


(3.20) 


Now  we  return  to  eq  (3.10)  and  note  that  R=a. 

.(Morse)  .  T (Morse)  6  aD  ,  ,n 

1 20  l22  =  - - —  exp[a(R0-a  )  J  1 


We  have 

-  exp(a(R0-a)l 

(3.21) 


Thus  , 


.  (1, axial)  _  r2  _al)  a(R0-a)(. 

*M2  =  2  2X  e  1 


ea(R0 -a)) 


(3.22) 


Th i s  quantity  is  greater  than  zero  only  as  long  as  R0  <  a.  For 
R0  >  a,  motion  at  the  centre  of  the  ring  is  not  stable.  The  mass 
cannot  remain  at  the  centre  of  the  ring  in  this  situation.  It  must 
move  to  a  position  of  equilibrium  elsewhere  in  the  system. 

For  the  transverse  motion  at  the  centre  of  the  ring,  we  find 

^(l.trans)  =  r2  2a?|)  e«(R.-u) 

(3.23) 

This  expression  depends  upon  two  dimensionless  parameters  which  may 
he  chosen  from  aRo ,  aa,  and  a/Ro. 

At  the  second  position  of  equilibrium  along  the  axis  of  the 
ring,  we  find  for  eq  (3.15) 


;u(Ro-a) (l.1/2rta) 


1 

7 


l 

RTa/ 
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(2, axial)  r2 
M2  7 


2ot2D(l  -  a2/Rg) 


(3.24) 


and  for  eq  (3.16) 


(2,trans)  r2 
M2  7 


(3.25) 


F.ach  of  these  motions  is  stable.  Indeed,  for  a2/Rj;  =  2/3,  the 
resultant  harmonic  potential  is  spherical .  This  particular  point 
lies  at  an  altitude  h  =  0.71a  above  the  plane  of  the  ring. 

For  the  Buckingham  exp-6  potential,  the  situation  is  more 
complicated.  Because  the  potential,  as  we  have  listed  it,  does  not 
depend  on  a  diatomic  equilibrium  separation  R0 ,  it  is  not  as  simple 
here  to  identify  the  condition  for  instability  at  the  centre  of  the 
ring  or  elsewhere. 

The  radial  equilibrium  condition  I i j (R  )  =  0  is  found  to  be 


-  R  R 

ABe  e  -  6C/R2  =  0. 


(3.26) 


The  general  second  order  term  is 


AB2 (1  - 2/BR) e 


-BR 


30C/R8  +  2P2 (cos0D)P2 (cos0  ) 

k  r 


x(AB2 (l-l/BR)e'BR  -  48C/R8) 


(3.27) 


Again,  for  the  consideration  of  equilibrium  at  thecentre  of 
the  ring,  R=a.  For  the  axial  displacement, 
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( 1  , axial )  v_z  1 

'B2  "27 


Here  we  find  that  if 


-oR 

AR  c  >  1 8  C  /  a 7 , 


6Sc-“»  .  i»JL  ' 

3  „  8 


(3.28) 


the  axial  disnlacement  at  the  centre  of  the  ring  will  he  unstable, 
For  the  transverse  displacement  ,  we  find 


( 1  ,  t  rans )  _  r2  1 


R2 


3 


iAR2  (3-5/aB)e'aR  -  S4C./a8 


(3.29) 


Vo r  the  second  position  of  equilibrium,  outside  the  ring, 
but  along  the  axis,  we  find 


R2 


r2  1  ad2  -B Rr 
j  ^AR  c  c 


0-7/BRe)  +  20-9/BRc)P2(cosO  )P2(cosPr  1 

e 


(3.301 


where  we  have  used  cq .  (3.26)  to  eliminate  (!/R^  terms.  };or  an 

axial  displacement,  we  find 


.12, a  x  i  a  1 
rR2 


r2  l.„2  -BRp 
=  T  *ARe  c 


3  - 


2  5 

im 


3  — (1-9/BR  ) 

R2  c 

e 


(3.31) 


For  the  transverse  displacement,  we  find 


(  2 , t  rans )  =  r?  lAg2 e ' BRe 


'^B 7 


r 


— (1+9/BR  )  -  2/PR 
2R2  c  c 


(3.32) 


In  these  cases,  it  is  necessary  to  test  numerically  for  equilibrium 
with  the  use  of  specific  values  of  the  quantities  A,  B,  and  C. 
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If  the  Coulomb  potential  is  added  for  a  charged  system,  it 
is  necessary  to  add  terms  of  the  form 


*c(R+r) 


q|Pjl(cos0R)P]l(coser) 


(3.33) 


in  which  r<(>)  is  the  lesser  (greater)  of  R  and  r.  The  addition  of 
terms  to  the  potential  clearly  makes  the  task  of  finding  stabilities 
more  difficult.  However,  the  formalism  does  allow  one  to  proceed 
directly  to  answer  such  questions  as  those  pertaining  to  stability. 


b.  Polygonal  distributions  of  discrete  sources 


We  write  the  general  expansion  of  an  arbitrary  scalar  function 
for  a  collection  of  sources  as 


C0 (R+r) 


/JW 


00  *  /s 

l  (rn/n!)J(-i)fc+nA  JP  (Rj.r)I  (Rj) 
n=0  St,  njij  x,  l  nx  i 


(3.34) 


where  the  summation  over  I  is  considered  over  sources.  For  a  regular 
polygon,  all  Rj  are  the  same,  R.  The  angular  distributions  with 
respect  to  the  centre  of  gravity  of  the  distribution,  of  course,  differ. 
Consider,  therefore,  the  sum  over  unit  vectors  R j : 


A  A 


l  pdRrr>  ’  HTT  l 

I  m=  -  £  l 


■  1  {feifr  P"(c°ser)Pm(coSeR)<=-i,',*r  | 


im<t>nT 

e  I 


(3.35) 


For  a  regular  polygon,^ 
discrete  sources,  and  n 


I 


=  2un/N,  where  N  is  the 
0,1,2 . N-l.  Thus  , 


total 


number  of 
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N-l 

I  exp  ( i  2nmn/N)  =  NS  ...  (3.3(>) 

n=0  ra»u 


Hence , 


t;0(r)  =  /4i  N  ^^(rn/n!)[(-i)£+nAnJ,P£(coseR)P^(cos0r)In£(R)  (3.57) 

where  each  0^  is  the  same  for  a  plane  of  sources. 

L-q  (5.37),  apart  from  the  number  of  discrete  sources  N,  is 
exactly  the  same  as  (3.3).  Hence,  all  conditions  of  stability  fournd 
for  the  continuous  ring  apply  equally  to  the  regular  polygonal 
distribution  of  sources. 

c.  Cages  of  sources:  regular  polyhedra 

We  now  consider  regular  polyhedral  cages  of  sources.  These 
cages  are  constructed  easily  by  distributing  sources  at  the  vertices 
of  tetrahedrons,  octahedrons,  etc.  The  tetrahedron,  for  example,  can 
be  realized  in  terms  of  a  base  triangle  with  a  fourth  source  located 
at  an  appropriate  position  on  the  axis  which  passes  through  the 
centre  of  the  triangle.  An  octahedron  is  realized  similarly.  The 
plane  of  sources  defines  a  square.  Two  additional  sources  arc  placed 
axially. 

in  order  to  investigate  the  stabilities  of  these  systems,  we 
make  use  of  the  expressions  for  for  polygons  together  with  as  many 
additional  terms  for  the  axially  located  sources  as  needed.  In  the 
following,  we  consider  the  tetrahedron  and  octahedron  as  the  most 
typically  encountered  examples.  More  complicated  systems  easily  are 
generated . 
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Considering  first  the  tetrahedron,  we  construct  a  base  triangle 
of  sources.  In  order  to  examine  the  mechanical  stability  of  a  test 
particle  located  at  the  origin  of  the  tetrahedron,  we  consider  a  point 
above  the  plane.  Thus, 

t2TCa)  =  ^  (- 1 2  o  CR)  +  2(l-3a2/2R2)P2(cos0r)I22(R)] .  (3.58) 

We  locate  the  fourth  source  on  the  z-axis  at  a  distance  R  from  the 
centre  of  the  tetrahedral  origin.  Thus, 

t,T(b)  =  /4¥  |2[-  2  o (R)  +  |P2(coser)I22(R)).  (3.39) 


The  sum  of  these  two  terms  gives 


tZl 


jl  2  o (R)  ♦  |(4-9a2/2R2)P2(cos0r)l22(R) 


(3.40) 


tf,  for  example,  we  consider  the  locations  of  the  sources  on  the 
triangle  to  be  (0 , 2/2R/3 , - R/ 3)  and  (+/271R, - /2R/3 , -R/3) ,  then  it  is 
easy  to  see  that  a2/R2  =  8/9.  Hence, 


t 


2T 


yl 2  o (R) 


(3.41) 


for  the  octahedron,  on  the  other  hand,  the  base  plane  of  (4) 
sources  contains  the  centrum  of  the  system.  Thus,  we  consider  the 
expression 


t 


(1) 

20(a) 


/41T 


r 4 
1 


r!2o(R)  "  jP2  (COSG^.)  I  2  2  (R) 


(3.42) 
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for  displacements  about  the  centre.  A  second  term  for  the  additional 
two  axial  sources  is 


t20(b) 


rl20  (R)  +  jP2  (COS0r)  I22  (R) 


The  sum  is 


(3.43) 


In  both  the  cases  of  the  tetrahedron  and  octahedron,  the  second 
order  Legendre  polynomial  disappears  as  is  required  by  arguments 
of  symmetry . 

It  is  useful  now  to  investigate  these  results  with  the  use  of 
specific  functional  forms.  The  use  of  the  Morse  potential  yields 
physically  transparent  and  interesting  conditions  for  stability. 

Tor  the  Morse  potential,  for  example,  we  can  write 


r 2  8 

^2M(T)  =  1  exP[<*(Ro-R)  ] 


2  exp [a (Ro - R) ] (1-1/aR)  -  1 


+  2/aR  . 


(5.45) 


If  R  =  Ro,  the  force  constant  for  an  encaged  particle  is 

kT  «  £a20,  (5.46) 

a  result  which  we  have  found  elsewhere.  Similarly,  the  force  constant 
for  the  displacement  of  a  test  particle  about  the  centre  of  an 
octahedron  of  sources  is 


-20- 


kQ  =  4  a 2  D .  13.47) 

In  both  of  these  cases,  it  is  clear  that  the  motions  are  mechanically 
stable . 

It  is  possible  to  consider  the  general  question  of  the  stability 
of  displacements  about  the  centre  of  gravity  and  symmetry  in  these 
structures.  From  eq  (3.19)  we  note  that  when 

2  exp[a(R0-R) ] (1-1/aR)  -  1  +  2/aR  <  0,  (3.48) 

a  displacement  is  mechanically  unstable.  The  condition  (3.48) --for 
the  Morse  potent ial -- reduces  easily  to  the  following  inequality: 

R  >  R0+  — £n2 .  (3.49) 

a 

When  this  inequality  is  satisfied,  instability  is  guaranteed. 

Such  conditions  of  instability  may  be  observed  in  actual  systems. 
Weaver  (personal  communication)  has  observed  that  the  attempted 
reduction  of  lithium  2.1.1  cryptate  to  form  the  encrypted  atom  is 
clectrochemically  irreversible.  The  gross  difference  in  size  of  the 
atom  as  compared  to  the  cation  suggests  the  system  may  be  unable 
to  accomodate  the  atomic  species.  If  the  encrypted  species  is 
mechanically  unstable  in  that  location,  then  some  form  of  tunnelling 
out  of  the  crypt  would  be  expected  rapidly  to  occur. 


-  21  - 

d.  A  truncated  Taylor  series  versus  the  exact  potential:  a  comparison 

We  conclude  this  section  by  comparing  a  truncated  Taylor 
series  for  the  Morse  potential  against  the  potential  itself.  We 
consider  again  the  tetrahedral  cage  of  sources,  each  of  which  is 
described  by  the  Morse  potential. 

The  choice  for  the  origin  of  the  Taylor  series  is  the  point  of 
equilibrium  at  the  centre  of  the  tetrahedron.  The  four  vectors 
(l<i<.4)  from  the  centres  of  force  to  that  point  have  the  same  length  R. 
Thus,  only  the  angular  factors  differ  in  the  four  individual  Taylor 
series.  The  superposition  of  the  Taylor  series  requires  the  considera¬ 
tion  of 


4  »  . 

I  VAR,  •  r )  (3.50) 

i  =  1  *  1 

in  the  manner  previously  discussed  [cf.,  the  last  subsection].  This 
resultant  angular  factor  can  be  evaluated  for  arbitrary  directions  of 
r.  This  angular  factor  (3.50)  vanishes  for  £  =  l  because  we  have 
chosen  the  point  of  equilibrium  to  be  the  origin.  It  also  vanishes 
for  £  =  2  because  the  system  has  no  quadrupole  symmetry.  We  were 
surprised  to  discover  that  it  also  vanishes  for  £  =  5,  although 
it  docs  not  vanish  for  any  other  value  of  £  up  to  12. 

In  order  to  demonstrate  the  non-spher ical  nature  of  the 
resultant  potential,  we  have  chosen  three  particular  directions  for 
r.  One  of  these  is  toward  a  vertex  of  the  tetrahedron,  i.e.,  directly 


toward  one  of  the  centres  of  force.  Another  is  toward  the  midpoint 
of  a  face  of  the  tetrahedron,  i.e.,  directly  away  from  one  of  the 


_  ?  o  _ 


centre's  of  force.  One  can  anticipate  that  these  should  he  respectively 

the  directions  of  the  most  rapid  and  least  rapid  increase  of  potential. 

The  third  direction  was  chosen  toward  the  midpoint  of  an  edge  of  the 
tetrahedron,  and  should  represent  any  intermediate  rate  of  increase 
of  the  potential.  The  exact  expression  for  the  resultant  potential 
can  also  be  readily  evaluated  along  these  directions. 

In  Figure  2  we  show  the  result  of  an  evaluation  for  the 

particular  choice  of  parameters  aR  =  4  and  aR0  =  (> .  We  note  the 

following  particular  features  which  are  revealed  in  the  figure.  (1) 

The  potential  is  spherical  and  harmonic  at  the  centre*.  (2)  The 
aspherical  and  anharmonic  contribution  is  well  described  by  n  =  3,  >=3 
out  to  about  r=0.1R  (one  third  of  the  way  to  the  face).  (3)  The 
resultant  potential  is  extremely  well  described  by  terms  up  to  n=o. 

Similar  results  are  obtained  for  an  octahedron  of  centres  of 
force.  in  this  case,  however,  the  resultant  potential  exhibits  an 
even  more  spherical  character.  The  exact  potential  is  effectively 
duplicated  with  the  use  of  terms  n=2,4,6.  For  the  tetrahedron,  in 
contrast,  n=2 , 3 , 4 , 5 , 6 ;  i.c.,  five  terms  arc  needed  to  give  a  good 
fit.  In  both  cases  the  actual  potential  functions  are  well  represented 
by  a  small  number  of  terms. 

4.  Summary  and  conclusions 


We  have  used  a  symmetry -adapted  form  of  the  Taylor  scries  in  order 
to  investigate  the  nature  of  the  stabilities  predicted  for  several 
mechanical  systems.  In  particular,  we  have  examined  rings,  polygons, 
and  polyhedra  all  of  which  share  similar  mathematical  properties. 
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The  examination  of  the  system  of  a  continuous  circular  distribu¬ 
tion  ol  matter  which  interacts  with  a  point  mass  showed  two  fundamental 
positions  of  equilibrium:  one  is  at  the  centre  of  the  ring  and  the  other 
is  on  the  axis.  We  established  the  conditions  for  the  points  of 
equilibrium  to  be  stable.  Next,  we  showed  that  a  regular  polygon  is 
mathematically  similar  to  the  continuous  ring.  The  conditions  for 
equilibrium  and  stability  are  the  same  for  both  systems.  Finally, 
we  showed  that  questions  of  stability  could  be  answered  for  regular 
polyhedra  by  making  use  of  the  results  for  the  polygons.  It  is 
possible  to  show,  for  the  cases  of  the  tetrahedron  and  octahedron,  that 
under  certain  conditions  even  an  enclosed  point  of  equilibrium  at  the 
centre  of  a  regular  solid  can  exhibit  a  mechanical  instability  with 
respect  to  displacements  about  that  point. 

We  believe  these  results  to  be  generally  useful  for  a  number  of 
applications.  The  point  of  view  which  we  adopt  in  carrying  out  our 
analyses  is  that  one  mass  be  examined  with  reference  to  the  static 
collection  of  surrounding  masses  (sources  of  force).  For  a  number 
of  systems,  it  is  sufficient  merely  to  focus  attention  on  a  single 
particle,  such  as  a  solvated  ion.  When  this  limited  view  is  too 
restrictive,  it  is  quite  simple  successively  to  look  at  each  mass 
in  a  collection  of  interacting  masses  and  develop  a  succession  of 
Taylor  series.  Thus,  in  this  manner  the  realistic  harmonic  vibrational 
force  field  can  be  developed  for  a  system  such  as  a  molecule.  The 
subsequent  development  is  a  picture  which  is  familiar  as  vibrational 
spectroscopy.  For  scalar  potentials,  the  use  of  our  form  of  the 
Taylor  series  enables  one  quickly  to  calculate  harmonic  and  higher 
order  contributions. 

The  analyses  of  mechanical  stability  for  polygons  and  polyhedra 
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are  important  in  the  cons iderat ion  of  several  transport  processes  , 

The  surface  of  a  metal  or  semiconductor,  for  example,  at  which 
electrochemical  oxidation  or  reduction  can  take  place  can  he  viewed 
accurately  as  an  extended  surface  network  of  interconnected  polygons. 
Thus,  the  question  of  the  stability  of  a  particle  which  interacts 
with  several  polygons  is  important  in  the  determination  of  the  process 
of  surface  migration.  Simple  systems,  such  as  these  extended  polygons 
as  surfaces  in  vacua,  can  be  rendered  more  complex  by  considering 
additional  distributions  of  molecules  about  the  surface.  The  question 
of  the  stability  of  a  single  particle  enmeshed  in  such  a  complex  is 
still  a  resolvable  question.  The  system  can  be  modelled.  Positions 
of  equilibrium  can  be  established  easily.  And,  finally,  stability 
against  displacement  for  a  particle  can  be  determined.  All  of  this 
can  be  done  for  a  wide  variety  of  model  potential  energy  functions. 

As  an  example,  in  a  separate  paper,  we  consider  the  relative  stabilities 
of  solvated  ions  and  atoms,  and  the  implications  these  stabilities 
have  on  the  deposition  of  a  metal  atom  at  an  electrode  surface.** 

This  work  was  supported  in  part  by  the  Office  of  Naval  Research, 
Arlington,  Virginia. 
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F  igure  captions  : 


Figure  1,  The  geometry  for  the  determination  of  the  equilibrium 
condition  near  a  ring  of  source  for  the  Morse  potential. 

figure  2.  The  potential  energy  as  a  function  of  displacement  along 
the  radii  in  a  tetrahedral  cage  of  sources  Tor  the  Morse  potential. 

As  indicated  in  the  inset,  the  radii  are  directed  toward  a  vertex 
(V)  ,  toward  a  face  (F)  ,  or  toward  an  edge  (F.)  •  The  maximum  displace* 
meat  considered  is  one-third  of  the  distance  to  a  vertex,  i.e., 
the  distance  to  a  face.  The  dashed  curve  in  each  labelled  pair  is 
produced  by  the  Taylor  series  to  order  r*,  while  the  solid  curve 
represents  (within  the  width  of  the  line)  both  the  exact  super- 
pot  ision  of  the  Morse  potent  ials  and  the  Taylor  series  to  order  r{' . 


